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Abstract: We consider the class of Rivlin Ericksen fluids of complexity 2 satisfying the Clausius- 
Duhem inequality. For strictly isolated fluids, we show that, if the Helmholtz free energy 
is not minimum in equilibrium, then the initial perturbation is (under some conditions) 
not only amplified, but also that it increases indefinitely. These results extend the analysis 
given by Fosdick Rajagopal in 1982 for the second grade fluids, which are a special case of 
Rivlin Ericksen fluids of complexity 2. 
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1 Introduction 


It is well-known that even if the Navier-Stokes equations are considered to be a good model 
for many fluids in the nature, there are many cases for which these equations can not apply. 
More and more fluids, existing in nature or created by industries (for example, some cristal 
liquids, bitums or oils), exhibit what called “non Newtonian” properties as the Weisenberg effect 
(climbing on a rod), the die-swell, a.s.o., properties observed and measured in practice. Fluids 
governed by the Navier-Stokes equations can not have such behaviors. Indeed, these properties 
are described by the viscometric functions T12, Ti: — T22 et To2 — T33, where T = (T;,5)1<i,j<3 
denotes the stress tensor. Consider a simple shear flow, i.e., the motion of a fluid between two 
walls, one fixed, the other one moving with the speed kz, (with k being the shear stress). In the 
Newtonian case, it is easily seen that T42 is linear in k, and moreover, Ty; — T22 = To — T33 = 0. 

In the non Newtonian case, it is observed that Tj is generally not linear in k, T31 — T22 #0 
and To2 — T33 Æ 0. So, for non Newtonian fluids, one has to introduce other constitutive laws 
than the law leading to the Navier-Stokes equations (modeling the Newtonian fluids) 


T = -pI + 2uD, 


where u is the velocity field, p is the pressure, p > 0 denotes the viscosity of the fluid, J is the 
unit tensor and D is the symmetrized rate of deformations tensor 


Vu + (Vut) 


D =2(L + L’) = 3 l 


with L = Vu. 
The non Newtonian theory is developed in several different directions. Let us mention some 
examples of non Newtonian fluids for which mathematical studies have been done: visco-elastic 


Received date: 2009-01-11. Biography: K. Boudra, Male, Doctor, Associate Professor. Research 
field: PDE. 


762 CHINESE JOURNAL OF ENGINEERING MATHEMATICS VOL. 26 


fluids (cf. for instance, Guillopé and Saut!!] and the references herein), Reiner-Rivlin fluids 
(cf. for instance, Duvaut and Lions!?), Ladyzhenkaya!!! Cioranescu!)) and the differential-type 
fluids (cf. for instance, Bernard'4], Cioranescu and Ouazar!®!, Cioranescu and Girault!®). 

The constitutive law for differential-type fluids (also called Rivlin-Ericksen fluids) was in- 
troduced by Rivlin!” (see also Truesdelll®!). In this theory it is supposed that the material 
is homogeneous and isotropic, and that the stress tensor depends only on the gradient of the 
velocity via the Rivlin-Ericksen tensors, defined reccurently as follows 


Ay =L+ Lt = 2D, can An re © Ant + LT An-1 =. An—1L. 


Then, the constitutive law for a differential fluid of complexity n is given (see Noll and Truesdelll®!) 
by 
T= —pI + F(A;, A2, para »An)s 


where F is an isotropic function, i.e., satisfying for any orthogonal matrix Q, 
QF(A, Aa, eee , An)QT = F(QA QT, QAQ7, Baad ,QAnQ™). 


The general form of an isotropic function F is given in Wang!!°l. For example, for a fluid 
of complexity two with density p > 0, 


F = pA, + a Ap + aA? + 03.43 + a4(A1A2 + A2A1) + 05(A?2 Ag + Ap A?) 
+ag(AZA; + A142) + a7 (AZAR + A3A}), 


where u = u(p, Ai, A2) and a; = a;(p, A1, A2) with 1 < i < 7, depend on the following func- 
tions, all isotropic: p, trA:, trA?, trA}, trA2, trA3, trA3, trA; Ae, trA? Ao, trA3Aj, trA?A2. 
The function F can also depend on the temperature 0. When F is a polynomial of degree n, 
the fluid is called of grade n. In this framework, the constitutive law of a fluid of grade 1 is 


T = —pI + pA, 
which in fact, is the definition of a Newtonian fluid. For a fluid of grade 2, one has 
T = —pI + pA; + a1 A2 + a2 A?. (1) 


Taking into consideration the dissipation Clausius-Duhem inequality, a thermodynamical study 
due to Dunn et Fosdick!!!], proved that the coefficients in (1) have to satisfy the following 
conditions 


u20, (2) 
a, +a2 =0. (3) 
Suppose moreover, that the Helmholtz free energy is minimum at the equilibrium, one also has 


Q1 > 0. (4) 
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Under these conditions, the fluid is stable. Dunn et Fosdick!!! and Fosdick et Rajagopal!!2] 
have also proved that if a, < 0 and a; + a2 = 0, then the motion of the fluid is asymptotically 
unbounded. Such a behavior has never been observed experimentally, it is why one may consider 
that the only admissible fluids of grade 2 could be those where (2)-(4) are satisfied. 

The positive sign of a; gives rise to many controversial discussions in the litterature (see 
for instance, Truesdell!§], Dunn and Rajagopall']). It is why, weaker restrictions than a; > 0 
have been introduced later on. In particular, Dunn!" studied the case where the coefficients 
in (1), depend on the temperature 


T = -pI + pa, Aı)Aı + an (6, A) Ag + a2(6, Ai) A?. (5) 


Under the restrictive hypothesis that the free energy is minimum at the equilibrium, he showed 
that 
a1 (0, A1) = a1 (0, tr A3), (6) 


and 
Zz 


fot ¿dE > 0, VzeE (Oe), €>0. (7) 
0 
This class of fluids, under appropriate hypotheses, has reasonable mechanical and dynamical 
stability properties. Conditions (6) et (7) are less restrictive than (4) in the constant case. 

Our aim here is to study the instability and the asymptotic behavior of the fluids defined 
by the law (5) and to weaken the restrictions on the coefficients a et a2 mentioned above. 
Following the ideas of Dunn and Fosdick!!! and of Fosdick and Rajagopal!], we do not require 
that the free energy of fluids defined by (5), be minimum at equilibrium (as is the case in 
Dunn!"4)). 

Theorem 2.3 below shows (under some additional hypotheses), that if is sufficiently large 
and if a2 € L3/2(Q), then any initial perturbation is not only amplified, but it also increases 
indefinitely. This result generalizes Theorem 1 of Fosdick et Rajagopal!!2], where a; and a2 
were constants. 

Proposition 2.5 describes the asymptotic behavior of the solution as t + +oo. We show 
that if p is sufficiently large, then 


lim inf (sup |Ai|(¢)) # 0, 


which is contrary to the usual Newtonian behavior. This result generalizes Theorem 2 from 
Fosdick et Rajagopal!!2!. 

Let now consider the cas where a1 + a2 = 0 on Q x [0, +00). Theorems 2.6 and 2.8 below 
extend to the non constant case of a1 and a2, the corresponding results from [11] and Theorem 
17 of [12]. They exhibit an anomaly of the type |{u(¢) 2 > cı N(0)exp At, for A € [0, 2u/c1], 
where N is defined by (11) and c; is given by (9) below. 


2 Setting of the problem and main results 


Let 9 c R be a bounded rigid connected domain with a boundary sufficiently smmoth. It 
is filled with a homogeneous incompressible fluid obeying to (5). We suppose that we are given 


764 
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an initial data at time t = 0, that on the boundary ôN one has an adherence condition and 


that the exterior field of forces derives from a potential. 


The flow of this fluid is described by the system 


divT + pf = poe, in Q x [0,+00), 


div u = 0, in Q x [0,+00), (8) 
u(x,t) =0, for c€ OQ, Vt>0, 
u(x, 0) = uo(z), for xz EQ, 


where p > 0 is the density of the fluid. é 


In the sequel we use the following notation: 

|A| = VtrAA? for any tensor A, tr stands for the trace of the tensor, 

Co is the Poincaré constant associated to Q; 

T, is the set of symmetric tensors while 72? the set of symmetric tensors with zero trace. 
To study the field of velocities u defined in problem (8), we will be only interested in regular 


motions, the temperature will not play any role since the field of temperatures is considered 


uniform in the domain Q and constant with respect to any t > 0. 


Throughout the paper we make the following assumptions: 

1) The density p is a positive constante and to simplify, we set p = 1, 

2) The viscosity yu is a strictly positive constant; 

3) The response functions a) and az are sufficiently smooth on Q x [0, +00) and such that 
aıl0, Ay) = ax (6, |4112), 
a2(0, A1) = a2(0, |41|?), (9) 
aı is bounded and |a;| < cı, where c; is a strictly positive constant. 

Our first result is the following: 


Proposition 2.1 For any fluid of complexity 2 defined by (5) and satisfying (8), one has 


|41]? 
d 1 
a f lul@dar + 2 J aı(0,£)dë + I nA de + | topti dza + <0) 
dt Q 2 Q Q 

0 


Let N(-): Rt + R be the function definded by 


|41]? 


Ne=-3 f ( J o1(0,E)d£)dz — fae. (11) 


Then we have the following result: 
Proposition 2.2 Let u be a velocity field satisfying (10). Suppose also that 


fawo < —2Co(1 +7)z, Vz>0, for some 7> 0. (12) 
0 
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Then, for every 2 and any field u(-,0), one has 


N(t)>0, Vt>0. (13) 
Suppose moreover, that 
u(-,0) #0, on Q, (14) 
and 
@2 is bounded in Q x [0, +00), (15) 
and also that 
lai +a2|>e on Ox [0, +00), for arbitrary € > 0. (16) 
Then 
(a) N(0)>0, (6) N(t)>0 on (0, +00). (17) 


In the same range of assertions, we have 
Theorem 2.3 Let u be a velocity field satisfying (10). Assume that (12), (14) and (16) 
are fulfilled. Moreover, if, 
a(t) E L972), Yt>0. 


Then, for a sufficiently large viscosity and for any domain Q, one has 
YM > 0, 3tm>0 such that f lai + af [Ay (F(a, tar)dar > M. 
Q 


Corollary 2.4 Let u bea velocity field satisfying (10). Assume that (15) and (16) hold. 
If the upper bound 6 of a, is strictly negative on Q x [0, +00), then in any arbitrary flow, for 
a sufficiently large viscosity and in a rather small domain Q, one has 


VM>0, Jtm>0 such that f |A1|3(z, tag)da > M. 
Q 


Proposition 2.5 Let again u be a velocity field satisfying (10). If (10), (14), (15) and 
(16) hold true, then, for a sufficiently large viscosity and for any Q, one has 


lim inf (sup 41) #0. 


Theorem 2.6 Let again u be a velocity field satisfying (10). Suppose that (12) and (14) 
are satisfied. If 
1 +2=0, on Q x [0,+00), (18) 


then 


lut)? > [(0) exp(At)Jcr, VA suchthat 0<\< 2h, 
1 


Corollary 2.7 Under hypotheses (14) and (18) and if the upper bound b of a; is strictly 
negative on Q x [0, +00), then 


rop = MOeeoo) VÀ such that o<a< t, 


Our last result is the next proposition. 


766 CHINESE JOURNAL OF ENGINEERING MATHEMATICS VOL. 26 


Theorem 2.8 Let us consider a fluid defined by (5) and system (8) with initial data 
satisfying (14). If Q c R?, then, independently of the choice of the sign of a; + ag (i.e., either 
lay + ag| > € or a1 + a2 = 0), one has 

2 
lut)? > [N(0)exp(At)Je:1, VA such that 0<A< E, 
1 


3 Proofs of the results 


3.1 Proof of Proposition 2.1 
It is sufficient to take the scalar product of the first equation in (8) with u, to integrate over 
Q and to use the diveregence theorem to obtain the following energy equation 


zg Pit fT Laz = f u: Tndo + | ufar. (19) 


Since f derives from a potential, with the help of (5), inequality (19) becomes 


jg [wrens fT- Ldr =0, Vt>0. (20) 


Making use of (5) and taking into account the symmetry of the stress tensor T’, one gets 


1 
T-L= 3T -A= —-pL+- 5 (uli? + Ial, Al 2) Ai? + (ay + az)trAĵ), (21) 


where we used the fact that 


Åz $ Ay = fap + trA?. 


z dt 
Substituting (21) in (20) and using again the fact that div u = 0, gives 


|41}? 


d 
T T E Í o0, e)d£+ | ulAiPar+ f (a1 + an)trAĝäe =0, 
0 


and this concludes the proof of Proposition 2.1. 
3.2 Proof of Proposition 2.2 
3.2.1 Proof of inequality (13) 
From (11) we can write 
|411? 


1 
N= =5 f J o1(0,£)dgde -Co | JAiPaz+Ca | lAde- f lupaz. (22) 
0 


[laPar=2 f |\Vul?da, 
Q Q 


since div u = 0 and u(-,t) = 0 on NQ. Applying the Poincaré inequality yields 


Observe now that 


|41]? 


1 2 
N() 2-5 f J on (6, €)déde — Co | |Aı|?dz. (23) 
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Due to hypothesis (12), one can estimate again the second term in the right hand side of (23) 
to get 


N(t) > 2(1 +nCn f |A;|?da ~ Co f |A;|?da, 
Q 2 


so that there is some ¢ > 0, such that 
N(t) > Ca f |A;|*da, Vt>0, 
Q 


and this ends the proof of (13). 

3.2.2 Proof of inequality (17)(a) 

Since u(x, 0) # 0 on Q, there exists £o € Q such that uo(xo,0) # 0. The continuity with 
respect to x of the velocity field u, implies the existence of a neighborhood Q% of £o, Qo C 2 of 
non zero measure, such that u(x,0) 4 0 for x € Qo. Consequently, there also exists a subset of 
Qo of non zero measure, on which A;(-,0) Æ 0 too. Indeed, if this is not true, as 2 is smooth 
enough and connected, the smoothness of OQ and the adherence condition of the fluid to the 
boundary of 2, would imply 


u(-,0) EC) and u=0 on Q, 
which is in contradiction with hypothesis (14). We deduce that we must have 
nca f |A,|?da >0, for t=0, 
Q 


when N (0) > 0. 
3.2.2 Proof of inequality (17)(b) 
The proof is based on the following result, due to Fosdick et Rajagopall’l. 
Lemma 3.1 For every A in Tf and for any a € R, one has 


-lA < atrA? < liar. 


Notice that by definition and thanks to the fact that div u = 0, one has A; € Tf. Now, let 
us consider the set 
E = {t € [0, +œ) | N(t) = 0}. 


If we suppose that it is not void, then it is easily seen that 


N(t) =0 <> Aj(z,t)=0 on Q. (24) 


f |Ai[2da = 2 f |Vuļ?dz,. 
2 2 


so that u(-,£) = 0 on Q for any t > 0, this being a consequence of the hypothesis that the 
motion is smooth, of the fact that Q is connected and of the boundary condition u = 0 on 89. 
Consequently, N (t) = 0 for any t > 0. 

On the other hand, if to is the lower bound of the set E, then N (to) = 0. It follows that 


Aj(-, to) =0 on Q. 
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As u = 0 on OQ, the smoothness of u and of 3N imply that A (-, t) belongs to the class C(Q). 
Consequently, A;(-,¢) admits a maximum on 2, maximum given by 


sup A;(z, ¢). 
Q 


Now, let cp be such that Ja; + a2| < c2 on Q x [0, +00). By using Lemma 3.1, one can 
estimate the following function 


4 vey = f ulAidz+ f (os + az)trAŞdz. (25) 
dt Q Q 
Indeed, 
d 2 
Aa > _ 2. 
ENO >u | MiPaz- Z f Afas, 
so that j 
C2 2 
SNO > [u - S supa] | Mafaz (26) 


But to is the smallest time for which A; (-, to) = 0 on Q (see (24)). As a consequence A,(-,t) #0 
on Q for any t such that 0 < t < to. It follows that for any t with 0 < t < to, there exists 
NQ: C 2 of non zero measure, on which A(x, t) Æ 0. 

Moreover, the continuity of Aı(-, -) allows us to find an interval [to — ô, to) such that 


0< sapje < wi, Vte [to — ô, to). (27) 


Recalling that 
A(z, t) #0, YTE UMN , 


t€[to—64,to) 
we get 
J |41|?dz >0. 
Q 
Therefore, by making use of (26) and (27), 
d 
J N(t) > 0, Vt E€ [to — 6, to). 


This implies that the function N is strictly incresing on [to — ô, to). Since N is strictly positive 
on [to — ô, to), one has 

lim N(t) > 0. 

ito 


Due to the continuity of N, one has N(tg) > 0, which is contradictory to N (to) = 0. So the set 
E is empty, and therefore N(t) is strictly positive for t > 0. 

3.3 Proof of Theorem 2.3 

Let À be an arbitrary real number. From (25), one has 


INe - AN(t) > ( (u- a) )f [Aida — J oaia] 3dr +d f lul@dx. (28) 


Choosing A such that 
2 
O<A< ae (29) 
cy 
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inequality (28) implies 
d 


TNE) -AN(t) > f fea eal Afan. (30) 


Suppose now by absurd, that 
4M>0 suchthat Vt>0, f lay + a2| |41| dz < M. (31) 
By virtue of hypothesis (16), 
[Faz < x, Vt>0. (32) 


Replacing (31) into inequality (30), yields 


O — \N(t) > -F 
Thus, 
£ (N (t) exp(—At)) > -5 exp(—At). 
Therefore 
NG) = (N0) - 5a) exp(At) + Tr 
ans sọ, 
NG) > (N0) - =a) exp(At). 


Since N(0) > 0, choosing u sufficiently large, according to (29), one can also take A sufficiently 
large, in order to get 


M 
N(0)- 57 >0. 
This implies an exponential growth of N(t) when t — +00. Then from definition (11) of N, we 
have 
|411? 
% J ai(9,f)d€dx > œ for t— +00, (33) 


0 
whence, by Holder inequality and since a; is bounded 


f |A;|?da = œ for t— +00. 
2 


But this result is in contradiction with (32). Therefore, (31) does not hold and this ends the 
proof. 

3.4 Proof of Corollary 2.4 

Choose Q in order to have 


Ca < 2(1 +7) 


for an arbitrary 17> 0. 


Hypothesis (12) is then satisfied and the corollary is immediate. 
3.5 Proof of Proposition 2.5 
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Proposition 2.2 states that N(0) > 0 on [0,+00). Then definition (11) leads to the following 


inequality 
|411? 


| a1(0,é)d€dz <0, Vt>0. 
a 
It follows that there exists Qo C Q of non zero measure, such that 


|411? 

J a1(0,é)dédz #0, Vt>0, 
Q 

0 


hence, 
|Ai|(z,t) #0 on Qo, Yt>0. 


Consequently, 
sup|A;|(¢) >20, Vt>0. (34) 
2 


We now prove that the strict positivity of 


sup |A;|(¢) 
Q 


is conserved as t — +00. Indeed, in view of inequality (34), only two cases can occur. 
Case 1 For any T € R*, there exists T > T such that 


sup |Ai|(r) 2 eae 


and so, 
lim inf sup |Ai|(t) 4 0. 
t—œœ gQ 


Case 2 There exists T > 0 such that for any t > T, one has 


0< sup | Ai|(t) < — n 


Then, by (27) et (29), one has 


|41]? 
iminf f | J aı(ð, aglae — f \u|?da > 0. 


Since a; is bounded, 
lim inf | |A,|?da > 0. 
t= Q 


On the other hand, as L°(Q) C L?(Q) with continuous imbedding, it follows that 
0< liming f |Ai|?dz < k liminf sup |A,|(¢), 
t—+00 Q t—++00 R 
where k is a positive constant. Consequently, 


lim inf sup |41 |(t) > 0, 
t—+00 Q 
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and the proof of Proposition 2.5 is now complete. 
3.6 Proof of Theorem 2.6 
From (25) and hypothesis (18), one has 


d 2 
SN) =H J Ai Pde, 


so that the function 


uf |Ay|?dx 
2 


is increasing on (0, +00). Let A be an arbitrary real number and introduce the following function 


F(t) = SNe) ~ \N(t). 


Then, recalling definition 11, it is immediately seen that 


|Aa|? 
F=f Pa- f | | an(0,)delde +A f luas. 


Since the function a, is bounded, if we choose À in the interval [0, 2), one has 
F(t)>0, Vt>0, 


and therefore, 
N(t) > N(O)exp(At), Yt>0. 


By Proposition 2.2, N(0) > 0, and so the function N has an exponential growth. 
From the expression of N, we get the following inequality 


a J |Ai?dz > N(0)exp(àt), Vt>0. 
Q 
But 
f {Ay |?da = 2f |Vul?dz, 
Q Q 


N(0) 


so that 


lull? > exp(At), Vt>0. 


3.7 Proof of Theorem 2.7 

The result is a simple application of Theorem 2.6. 
3.8 Proof of Theorem 2.8 

Applying Hamilton-Cayley theorem, we get 


A? + IdetA = 0. 


Multiplying by A and taking into account that trA = 0, it follows that trA? = 0. Then, the 
theorem follows from all the preceding results. 
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